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1. Introduction 

Low dimensional strongly coupled electron-phonon systems like MX-chains, ferroelec- 
tric perovskites, conjugated polymers, or organic charge transfer salts exhibit a remark- 
able wide range of strengths of competing forces and, as a result, physical properties [1]. 

Most notably quasi one-dimensional (ID) materials are very susceptible to structural 
distortions driven by the electron-phonon (EP) interaction. Probably the most famous 
one is the Peierls instability [2] of ID metals: As the temperature is lowered the system 
creates a periodic variation in the carrier density by shifting the ions from their symmet- 
ric positions. For the half-filled band case this so-called charge density wave (CDW) is 
commensurate with the lattice, the unit cell doubles, and the system possesses a spon- 
taneous broken-symmetry ground state. Since a static dimerisation of the lattice opens 
a gap at the Fermi surface the metal gives way to a Peierls insulator (PI) [see Fig. 1]. 

The on-site Coulomb interaction, on the other hand, tends to immobilise the charge 
carriers as well by establishing a Mott insulating ground state. The Mott insulator (MI) 
exhibits strong spin density wave (SDW) correlations but has continuous symmetry and 
therefore shows no long-range order in ID. Then, of course, the question arises, whether 
the PI and MI phases are separated by one (or more than one) quantum critical point (s) 
at T = 0, and if so, how the cross-over is modified by quantum phonon effects. 
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Fig. 1. - Peierls scenario: A gap 2A opens in the electronic band structure E(k) [density of 
states p{E)] of an ID metal if, as a result of the EP coupling, a static lattice distortion occurs, 
implicating a new lattice period 2a in real space. 



The challenge of understanding such quantum phase transitions has stimulated intense 
work on generic microscopic models of interacting electrons and phonons like the Holstein 
Hubbard model (HHM): 

(!) h = -t c t<r c 3° - 9^° J2( b l + b ^ n ™ + w ° h \ h i + U Y1 n ^ n n ■ 
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Here c\ a {c ia ) denote fermionic creation (annihilation) operators of electrons with spin 
a =|, J. on a ID lattice with N sites, ni a = c| c^, and b\ (b^) are the corresponding 
bosonic operators for dispersionless optical phonons. 

The physics of the HHM is governed by three competing effects: The itinerancy of 
the electrons (oc t), their on-site Coulomb repulsion (oc U), and the local EP coupling 
(oc g). Since the EP interaction is retarded, the phonon frequency (loq) defines a further 
relevant energy scale. This advises us to introduce besides the adiabaticity ratio (uio/t) 
two dimensionless coupling constants (u — U/At and A = 2e p /2t or g 2 = e p /u ). 

In the single-electron case, the Holstein model [3] has been studied extensively as a 
paradigmatic model for polaron formation (see, e.g., Ref. [4]). Here e p gives the polaron 




Fig. 2. - Schematic phase diagram of the ID Holstein Hubbard model. At half-filling, Peierls 
(left) or Mott (right) insulating phases may be favoured over the metallic state. In the case 
of localised electrons interacting via antiferromagnetic exchange and magneto-elastic couplings 
even a spin-Peierls distorted state can emerge (right, lower panel). 
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binding energy. The Hubbard model [5] , originally designed to describe ferromagnetism 
of transition metals, more recently has been in use in the context of metal-insulator 
transition and high-temperature superconductivity as the probably most simple model 
to account for strong Coulomb correlation effects. As yet there exits almost no exact 
(analytical) results for the full HHM (1). At least at half-filling, however, it has become 
generally accepted that the interplay of charge, spin and lattice degrees of freedom gives 
rise to the phase diagram sketched in Fig. 2. This scenario is supported by dynamical 
mean field investigations of the HMM, which become reliable at least in infinite spatial 
dimension [6]. Besides the properties of the ground state, the nature of the physical 
excitations is puzzling as well, especially in ID. While one expects "normal" electron- 
hole pair excitations in the PI phase (U = 0), charge (spin) excitations are known to 
be massive (gapless) in the MI state of the Hubbard model (A = 0). Thus, varying the 
control parameter u/X, a cross-over from standard quasi-particle behaviour to spin-charge 
separation might be observed in the more general ID HHM. 

The aim of this contribution is to establish this physical picture and the anticipated 
phase diagram of the ID HHM. For these purposes we employ the Lanczos exact diagonal- 
isation (ED) [7,8], kernel polynomial (KPM) [9,10] and density-matrix renormalisation 
group (DMRG) [11] methods (see Jcckclmann/Fchskc [12]). These numerical techniques 
allow us to obtain unbiased results for all interaction strengths with the full quantum 
dynamics of phonons taken into account. 

2. Luttinger-liquid Peierls-insulator transition 

2T. Holstein model of spinless fermions. - In a first step, let us neglect the spin degrees 
of freedom. The resulting ID spinless Holstein model, H = — i 2Z^<-i j> c I c j — fl^o + 
bi) n i + w oSt^i^»> exhibits a quantum phase transition from a metallic to an insulat- 
ing phase at half-filling (A = A/2) [13, 14], where the critical coupling X c (uio) — > for 
luq — > 0. In the anti-adiabatic (uq — » oo) strong EP coupling regime, the model can be 
transformed to the exactly solvable XXZ (small polaron) model [13], which shows a tran- 
sition of Kostcrlitz-Thouless type. Various variational [15], renormalisation group [16], 
world-line quantum Monte Carlo [13] or Green's-function Monte Carlo [17] methods were 
used to determine the phase boundary, within which significant discrepancies occur in the 
adiabatic intermediate coupling regime. More precise ED [18] and DMRG [19] techniques 
yields the phase diagram presented in Sec. 2'3. 

2'2. Luttinger-liquid parameters and charge structure factors . - Before we discuss the 
metal insulator transition in the framework of the Holstein model we will characterise 
the metallic and insulating phases in themselves. According to Haldane's Luttingcr liq- 
uid (LL) conjecture [20], an ID gapless (metallic) system of interacting fermions should 
belong to the Tomonaga-Luttinger universality class. Since the Holstein model of spin- 
less fermions is expected to be gapless at weak couplings g, we try to determine the 
(non-universal) LL parameters, K p (correlation exponent) and u p (charge velocity), by 
performing a large-scale DMRG finite-size scaling analysis. To leading order, the ground- 
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Fig. 3. - DMRG finite-size scaling of the charge gap A C (N) and the ground-state energy Eq(N). 
ED data included for comparison. The Table gives the LL parameters extracted from the scaling 
relations. Note that K p > 1 (K p < 1) in the adiabatic (anti-adiabatic) regime. 



state energy and the charge excitation gap of a finite system with TV sites scales as: 

(2) IT =£ » (00) -3T^' 

(3) Ac(A0 = E±(N) E (N) = ^±.L. 

Here £o(°°) denotes the energy density of the infinite system with N/2 electrons, and 
Eq(N) are the ground-state energies with ±1 fcrmions away from half-filling. 

The LL scaling relations (2) and (3) were derived for the pure electronic spinless 
fermion model only [21]. Figure 3 demonstrates, exemplarily for the adiabatic regime, 
that they also hold for the case that a finite EP coupling is included. The resulting LL 
parameters are specified in the Table. Interestingly the LL phase splits in two different 
regimes: For low phonon frequencies the effective fermion-fermion interaction is attractive 
(K p > 0), while it is repulsive (K p < 0) for high frequencies. In the latter region the 
kinetic energy (oc u p ) is strongly reduced and the charge carriers behave like (small) 
polarons. In between, there has to be a point where the LL is made up of (almost) non- 
interacting particles (K p = 1). The LL scaling breaks down just at the critical coupling 
g c (u)o/t), signalling the transition to the CDW state. We found g^((Jo/t = 0.1) ~ 7.84 
and g%{u /t = 10) ~ 4.41. 

Figure 4 proves the existence of CDW long-range order above g c . Here the staggered 
charge structure factor 

(4) S c (tt) = i_£(-l);((n; - \){n i+j \)) 

i,3 



unambiguously scales to a finite value in the thermodynamic limit (N — ► oo). In contrast 
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we have S c (n) — ► in the metallic regime (g < g c ). Note that such a finite-size scaling, 
including dynamical phonons, is definitely out of the range of any ED calculation. 

2'3. Phase diagram of the Holstein model of spinless fermions. - Figure 5 gives the 
position of the LL-CDW phase boundary in the g-0J$ l plane. For weak EP interactions 
the system is a metal with LL parameters depending on both the coupling strength and 
the phonon frequency. Increasing uj at fixed g (which nonetheless means that the EP 
coupling parameter A becomes larger) the cross-over from an attractive LL (adiabatic 
regime) to a repulsive LL (anti- adiabatic regime) takes place at about uj^/t c± 1. Keep- 
ing loo fixed we enter the CDW phase at a critical EP coupling g c (uj )- ED, cluster 
perturbation theory [22] and projector-based renormalisation methods [23] reveals the 
softening of the optical phonon at the Brillouin-zone boundary, at least in the adiabatic 
regime, which can be understood as precursor effect of the gap formation. Note that 
in D = oo, the opening of the electronic gap is accompanied by the appearance of a 
low-energy phonon peak in the total phononic spectral function [24] . 

The CDW for strong EP coupling is connected to a Peierls distortion of the lattice, and 
can be classified as traditional band insulator and polaronic superlattice in the strong- 
EP coupling adiabatic (ujo/t <C 1) and anti-adiabatic (u>o/t 3> 1) regimes, respectively. 
Extremely valuable information about the CDW state can be obtained by analysing the 
regular part of the optical conductivity, 

(5) o-M = £ l ^ |itEj( t J+1 V 1+lC,)l ^ ) ' 2 *> - - ^ ' 

E m -E a 

which is connected due to finite-frequency optical transitions to excited quasi-particle 
states \^ m )- In Eq. (5), cr rcg (w) is given in units of we 2 and we have omitted an 1/AT 
prcfactor. The evaluation of dynamical correlation functions like a rc& (uj) can be carried 
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Fig. 5. - Ground-state phase diagram of the ID half-filled Holstein model of spinless fermions, 
showing the boundary between the LL and CDW states obtained by ED and DMRG approaches. 
The dashed line gives the asymptotic result for the XXZ model. Left insets show the LL 
parameters u p and K p as a function of the EP coupling g in the metallic regime; right insets 
display for a six-site chain the regular part of the optical conductivity a rcs (uj) (dotted lines) and 
the integrated spectral weight S rce (uj) = J u dio'a rcs (uj') (solid lines) in the CDW region. 



out by means of the very efficient and numerically stable ED-KPM algorithm [10]. The 
optical absorption spectra shown in Fig. 5 elucidate the different nature of the CDW 
for small and large adiabaticity ratios. In the adiabatic region the most striking feature 
is the sharp absorption threshold and large spectral weight contained in the incoherent 
part of optical conductivity. In the anti-adiabatic regime the CDW is basically a state 
of alternate self-trapped polarons, which means that the electrons are heavily dressed by 
phonons. Since the renormalised band dispersion is extremely narrow, finite-size gaps 
are reduced as well. Therefore, A opt read off from Fig. 5 correctly gives the CDW gap. 

3. Peierls-insulator Mott-insulator transition 

31. Ground-state properties. - Now we include the spin degrees of freedom and ask 
for the effect of a finite Coulomb interaction. The ground state of the pure Holstein 
model (U = 0) is a Peierls distorted state with staggered charge order, i.e. alternating 
empty and doubly occupied sites, for g > g c (uJo) [25,26]. As in the Holstein model of 
spinless fermions, quantum phonon fluctuations destroy the Peierls state for g < g c . 
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Fig. 6. - Spin and charge struc- 
ture factors at q — ir in the half- 
filled ID 8-site HHM (1) with 
PBC for different u at A = 1. 
Squares denote ED results, filled 
circles show DMRG calculations 
with m = 600 target states and 
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The charge structure factor, S c (n) [cf. Eq. (4], and spin structure factor, 
(6) S s (n) = ± J2(-l) j (S*S* +j ) with S? = \{n t] - n t[ ) , 

shown in Fig. 6 for the full HHM, indicate pronounced CDW and weak SDW correla- 
tions provided u/X < 1. Increasing the Hubbard interaction u at fixed EP coupling A 
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Fig. 7. - DMRG study of staggered charge 
and spin (inset) structure factors as a func- 
tion of inverse lattice size in the PI and MI 
regimes, as well as at the quantum criti- 
cal point (QCP). Generally, m = 1000 and 
5 boson pseudo-sites were chosen for the 
DMRG algorithm, with the exception of 
the JV = 128 case where m = 800 had to 
be used for practical reasons with no qual- 
itative loss of accuracy. 
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Fig. 8. - Phonon distribution, \c™\ 2 , 
for typical model parameters char- 
acterising the MI, band insulator 
(BI) and bipolaronic insulator(BPI) 
Peierls phases (inset). Note that we 
have separated the Q — centre of 
mass phonon. 



and frequency ujq, the CDW correlations become strongly suppressed, whereas the spin 
structure factor at q = tt is enhanced. 

The results depicted in Fig. 6 are obtained for a rather small eight-site chain with 
PBC. In order to conclude about a possible existence of charge and/or spin long-range 
order we have calculated S c (ir) and S s (ir) for different system sizes and performed a 
finite-size scaling of our DMRG results (see Fig. 7). In the PI phase, S c (n) shows almost 
no dependence on the size of the system, indicating true CDW long-range order, whereas 
S s (ir) obviously scales to zero as N — ► oo. By contrast, in the MI regime our data 
provides strong evidence for vanishing charge but also spin order in the thermodynamic 
limit. Clearly the MI is characterised by short-ranged antiferromagnetic spin correlations 
but nevertheless the staggered spin-spin correlation function shows a slow (algebraic) 
decay at large distances. 

Figure 8 gives the phonon distribution function, i.e. the weights of the m-phonon 
state for different ground states of the HHM. First of all the results demonstrate that 
our Hilbcrt space optimisation and truncation procedure is well-controlled in the sense 
that states with larger number of phonons, as accounted for in the calculations, have 
negligible spectral weight. Of course, the number of phonons which have to be taken 
into account depends on the physical situation. Whereas the ground state of the MI is 
basically a zero-phonon state, multi-phonon states become increasingly important if u/X 
is reduced, i.e. in the PI state. 

3'2. Optical response. - The evolution of the optical conductivity cr rcg (w) going from 
the PI to the MI phase by increasing u/X is illustrated in Fig. 9. In the PI regime 
the electronic excitations are gapped due to the CDW formation. Excitonic gap states 
may occur in the process of structural relaxation. Because the PI ground state is a 
multi-phonon state, we find a gradual rise of the integrated spectral weight of a rcg (uj). 

At the QCP the optical gap A opt closes. Since we are in the non-adiabatic regime 
(luq ~ t), the lowest optical excitations have mainly pure electronic character in the 
vicinity of (u/X) c , i.e. the gap is closed by a state having large electronic spectral weight. 
Due to the selection rules for optical transitions the PI-MI transition necessarily implies 
a ground-state level crossing with a site-parity change. The site inversion symmetry 
operator P is defined by Pc\ a P^ — c ! N _ ia (N = An). We have explicitly verified that the 
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Fig. 9. - Optical conductivity in the Holstein Hubbard model (ED, N = 8, PBC). Dashed 
lines give the integrated spectral weight normalised by S reg (oo). S TCg (ui) / S TBg (oo) is a natural 
measure for the relative weight of the different optical absorption processes. Note that the 
critical coupling is identical to the point where S c (tt) sharply drops (cf. Fig. 6). 



parity is P = +1 (P = -1) in the PI (MI) phase. 

In the MI the optical gap is by its nature a correlation gap. The lower panel in Fig. 9 
shows clearly that CT reg (w) is dominated by excitations which can be related to those 
of the pure Hubbard model. In addition, phononic sidebands with low spectral weight 
appear. 

3'3. Photoemission spectra. - Next we study the spectral density of single-particle ex- 
citations associated with the injection of a spin-cr electron with wave number K, A~^ a (uj) 
(inverse photoemission (IPE)), and the corresponding quantity for the emission of an 
electron, A^ ct (w) (photoemission (PE)), where 

(7) A±» = X>ml4<>o>| a 8[w T {Et - E )] 

m 

with c\ ta = c' Ka and c Ka = c Ka . |-0 O ) , E refer to the ground state of the system 
with N c = N electrons and iV'm) (^m) are eigenstates (energies) of the (N c ± l)-particle 
system. Adding the spectral densities of (photo-) emission and absorption we obtain 
the spectral function = A^- a (uj) + A^^iuj), which obeys various sum rules and 

allows for a connection to angle-resolved photoemission spectroscopy (ARPES). 

Figure 10 displays the IPE and PE spectra for the HHM at the allowed wave numbers 
of our finite system: K = 0, ±7r/4, ±7t/2, ±37t/4, and it. To reliably monitor a band 
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splitting induced by the Hubbard and EP couplings it is necessary to guarantee that the 
Fermi momenta Kp = ±ir/2 are occupied, which is the case for N = 4/ (I integer, PBC). 

The most prominent feature we observe in the PI regime is the opening of a gap at 
K = ±tt/2. For the BI a rather broad (I)PE signatures appear. Within these excitation 
bands the spectral weight is almost uniformly distributed, which is a clear signature of 
the multi-phonon absorption and emission processes that accompany every single-particle 
excitations in the PI. The lincshapc reflects the (Poisson-like) distribution of the phonons 
in the ground state. The lower and upper band closely follow a (slightly renormalised) 
cosine dispersion. The situation changes radically if the insulating behaviour is associated 
with localised bipolarons forming a CDW state (see Fig. 10, upper right panel). Due to 
strong polaronic effects an almost flat band dispersion results with exponentially small 
(electronic) quasi-particle weight. Now the dominant peaks in the incoherent part of the 
(I)PE spectra are related to multiples of the (large) bare phonon frequency. 

If we increase the Hubbard interaction at fixed EP coupling strength the Peierls gap 
weakens and finally closes at about (u/X) c ~ 1, which marks the PI-MI cross-over. This is 
the situation shown in the lower left panel. Approaching the QCP, the ground state and 
the first excited state become degenerate. The QCP is characterised by gapless charge 
excitations at the Fermi momenta but perhaps should not be considered as metallic 
because the Drude weight in the case of a degenerate ground state is ill-defined [27]. 

If the Hubbard interaction is further increased, i.e. the Coulomb repulsion overcomes 
the attractive on-site EP coupling the electronic band structure becomes gapped again 
and a MI state develops (see lower right panel). The Mott-Hubbard correlation gap 
almost coincides with the optical gap A opt determined by evaluating the regular part 
of the optical conductivity for the same parameters. The form of the spectra, however, 
is quite different from PI case. Contrary to the BI phase in the MI regime the lowest 
peak in each K sector is clearly the dominant one. Then the dispersion of the lower 
(upper) Hubbard band can be derived tracing the uppermost (lowest) excitations in each 
K sector. Due to the finiteness of our system and the rather moderate value u = 1.5, PE 
(IPE) excitations with K = ±37r/4 and ir (K = ±7r/4 and 0) have still finite spectral 
weight. Since the spectral weight of the PE excitations with K > tt/2 is expected to 
vanish as N goes to infinity for u 3> 1, the lower Hubbard band will be completely filled, 
and consequently the system behaves as an insulator at T = 0. As a result of the coupling 
to the phonon system the electronic levels in each K sector split, creating phonon side 
bands. The distinct peaks are separated by multiples of the bare phonon frequency and 
can be assigned to relaxation processes of the Q — phonon modes [23] . The number 
of phonons involved is controlled by g 2 . S^ a {uj) shows clearly that the total spectral 
weight of the resulting excitation bands equals the weight of the respective electronic 
excitations in the pure Hubbard model. Interestingly, mediated by Q ^ phonons, 
there appear "shadows" of the bands composed of a dominant electronic excitation and 
phonon satellites in a certain K sector (e.g. K = ±7r/2) in the other K sectors (e.g. 
K = ±7r/4). These signatures give rise to a weak "breather- like" excitation being almost 
dispersionlcss in the Brillouin zone. The formation of quantum breathers was proposed 
by W. Z. Wang et al. [29], but without doubt has not been detected so far. 
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Fig. 10. - Wave-number-resolved spectral densities for photoemission {A~^ a (uj); dashed (red) 
lines) and inverse photoemission (Aj^ CT (w); solid (black) lines) obtained by applying our ED- 
KPM scheme to the HHM. Shown are typical results obtained for the case of a Peierls band 
insulator (BI) at wo/( < 1 (left) and bipolaronic insulator (BPI) at uio/t S> 1 (right). The 
corresponding integrated densities S^ a (oj) are given by dashed lines. The lower panels give the 
PE and IPE spectra near the PI-MI transition point (u/ ~ 1) and in the Mott insulating state 
(it/A S> 1). Data for the pure Hubbard model (blue and green lines) were shifted by —(e p N e /N) 
and included for comparison. Results are taken from Ref. [28]. 



3'4. Many-body excitation gaps. - Since many-body gaps to excited states form the 
basis for making contact with experimentally measurable excitation gaps and can also 
be used to characterise different phases of the HHM, we finally determine the charge and 
spin gaps, 



(8) 
(9) 



A c = E+(l/2)+Eo(-l/2) - 2E (Q) 
A s = E (l) - E (0) , 
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in the limit of large u > 1 [30] . 



using DMRG. Here E^ ±} (S Z ) is the ground-state energy of the HHM at half-filling (with 
N c = N ± 1) particles in the sector with total spin-z component S z . 

The scaling of A c and A s is shown in fig. 11. Obviously, A c and A s are finite in the 
PI and will converge further to the same value as N — > oo. Both gaps seem to vanish 
at the QCP of the HHM with finite-frequency phonons, but in the critical region the 
finite-size scaling is extremely delicate. In the MI we found a finite charge excitation 
gap, which in the limit u/A>l scales to the optical gap of the Hubbard model, whereas 
the extrapolated spin gap remains zero. This can be taken as a clear indication for spin 
charge separation. 

4. Summary 

In this report we have addressed the important problem of quantum phase transi- 
tions in one-dimensional strongly coupled electron-phonon systems. As a generic model 
we analysed the Holstein Hubbard model at half-filling. Applying numerical diagonal- 
isation methods we obtained, by the use of present-day leading-edge supercomputers, 
basically exact results for both ground-state and spectral properties in the overall region 
of electron-electron/clcctron-phonon coupling strengths and phonon frequencies. 

For the spinless Holstein model we found that for weak EP couplings the system 
resides in a metallic (gapless) phase described by two non-universal Luttinger-liquid pa- 
rameters. The renormalised charge velocity and the correlation exponent are obtained 
by DMRG from finite-size scaling relations, fulfilled with great accuracy. The Luttinger 
liquid phase splits in an attractive and repulsive regime at low and high phonon fre- 
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quencies, respectively. Here the polaronic metal, realised for repulsive interactions, is 
characterised by a strongly reduced mobility of the charge carriers. Increasing the EP 
interaction, a cross-over between Luttinger-liquid and charge-density-wave behaviour is 
found. The transition to the CDW state is accompanied by significant changes in the 
optical response of the system. Most notably seems to be the substantial spectral weight 
transfer from the Drude to the regular part of the optical conductivity, indicating the 
increasing importance of inelastic scattering processes in the CDW (PI) regime. 

For the much more involved Holstein Hubbard model, with respect to the metal the 
electron-electron interaction favours a Mott insulating state whereas the EP coupling is 
responsible for the Peierls insulator to occur (see Fig. 12). True long-range (CDW) order 
is established in the PI phase only. The PI typifies a band insulator in the adiabatic weak- 
to-intermediate coupling range or a bipolaronic insulator for non-to-antiadiabatic strong- 
coupling. Our results for the single-particle spectra indicate that while polaronic features 
emerge only at strong EP couplings, pronounced phonon signatures, such as multi-phonon 
bound states inside the CDW gap, can be found in the Mott insulating regime as well. 
This might be of great importance for interpreting photoemission experiments of low- 
dimensional materials such as MX-chain compounds. The optical conductivity shows 
different absorption features in the MI and PI as well and signals that the quantum 
phase transition between these phases is connected to a change in the ground-state site- 
parity eigenvalue (of a finite HHM systems with PBC) . From our conductivity data we 
found evidence for only one critical point separating Peierls and Mott insulating phases 
in the Holstein Hubbard model with dynamical phonons. This differs from the results 
obtained in the adiabatic limit (o>o = 0), where two successive transitions have been 
detected for weak couplings u, A « 1 [26]. The Peierls to Mott transition scenario is 
corroborated by the behaviour of the spin- and charge excitation gaps. From a DMRG 
finite-size scaling we found A c = A s and A c > A s = in the PI and MI, respectively. 
The emerging physical picture can be summarised by the phase diagram shown in Fig. 12. 
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Fig. 12. - Sketch of the PI-MI quantum phase transition in the HHM. 
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